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000000 
000000 
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6-SEP=-19 MTHRTL.SRCJOTSPOWCDC MAR; 1 (1) 
1 -TITLE OTSSPOWCDCD =- D SOMPLERSIC. ** D COMPLEX*16 routin 
; eIDENT /1-002/ 3; File: OTSPOWCDC.MAR Edit: $8L1002 
¢ s Aue enecoqanenoseeereonsoncessececqnqqsnoneccesqqoensencseeesotseenecetersees 
**® w i 
0 5 [* COPYRIGHT (c) 1978, 1980, 1982, 1984 BY * | 
0 8 ;* DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. * 
33 18 :* ALL RIGHTS RESERVED. : 
$00 11 ie THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND copies * 
44 \¢ 3* ONLY IN ACCORDANCE WITH T OF suUC ICENSE D WITH * 
0 15 ;* INCLUSIO HE ABOVE COPYRIGHT ICE. THI ARE OR A OTHE * 
0000 14 ;* COPIES THEREOF MAY NOT BE VIDED OR OTHERWISE MADE AVAILABLE TO ANY * 
0000 15 ;* OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
4464 1 3* TRANSFERRED. 5 
3* ® 
0000 18 ;* THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 
0000 19 ;* AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT * | 
443 $9 i* CORPORATION. : 
0000 3 * DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS * 
0000 3* SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. * 
0000 «98 is : 
0000 6 MISIIIIIIIIIII IIIT TIT TTI I IIT TTTiTiTi ti tlii iti itt titi iTiiiiiti titi iil 
0000 7: 
ee 
444 $0 > FACILITY: Language support Library - user callable 
9900 3¢ > ABSTRACT: 
9000 34 : D COMPLEX*16 base to D COMPLEX®*16 power giving D COMPLEX*16 result. 
0000 36 : ENVIRONMENT: User Mode, AST Reentrant 
ane 
0900 33 : AUTHOR: Steven B. Lionel, CREATION DATE: 20-July-1979 
9000 41 : MODIFIED BY: 
0000 rk : 1-001 - Original. Adapted from OTS$POWCC version 1-003. SBL 20-Jul-1979 
0000 44 : 1-002 - Use general mode addressing. SBL 30-Nov-1981 | 


] 
tc 9 
TSSPOWCDCD - D COMPLEX*16 ** D COMPLEX*16 routine 16-SEP-1984 AX /VMS o V04-00 P 
9805 DECLARATIONS ‘ on see Oke oN 33: 3 YATURTE L.S SOTSPOWCDC.MAR: 1 veg 3 
0 4g .SBTTL DECLARATIONS | 
48 : INCLUDE FILES: | 
49 ; 
0 } 
13 
9 § ; EXTERNAL DECLARATIONS: | 
as DSABL GBL | 
5 “EXTRN RTs coExP ; Complex exponentiation 
0 § -EXTRN MTHSCOLOG : Complex Logarithm 
$8 ¥ SEXTRN OTSSMULCD.. R3 ; Complex multiplication 
00 9 : MACROS: 
00 60 ; 
4 61 
00 6¢ 3 
0000 635 ; EQUATED SYMBOLS: 
44 64 ; 
00000004 000 65 base = 4 ; base input - by value 
00000014 44 6 exp = 20 3 exponent input - by value 
$009 68 : OWN STORAGE: 
0000 69 ; 
0000 70 
0000 71 ; 
B00 i 3 PSECT DECLARATIONS: 
sai ++ 7% -PSECT _OTSSCODE PIC, USR, CON, REL, LCL, SHR, - 
000 75 EXE, RD, NOWRT. LONG 
0000 76 
} 


Dd 9 
TSSPOWCOCD = D COMPLEX*16 ** D COMPLEX*16 routine 16-SEP-1984 01:55:2 AX/VMS Macro V04-00 Pp 
#8 OTSSPOUCDED RS - D COMPLEX*16 ** D COMPL met 9133 ‘23 EMTHRTL. SREJOTSPOWCDC.MAR: 1 - ad | 
8 «SBTTL OTSSPOWCDCD_R3 = D COMPLEX*16 ** D COMPLEX®16 
9 ; FUNCTIONAL DESCRIPTION: 
§ ; OTSSPOWCDOCD_R3 evaluates the result of taking a complex base 
3 to a complex power. The ANS FORTRAN X3.9-1978 standard defines 
° ¢ : complex exponentiation as: 
0 6 : x # y = CEXP( y * CLOG(x)) 
3 8 3 where x and y are type D COMPLEX*16. | 
09 $' ; The arguments of OTSSPOWCDCD_R3 are CALL BY VALUE. | 
$8 3 > CALLING SEQUENCE: 
$4 3¢ ; power.wdc.v = OTSSPOWCDCD_R3 (base.rdc.v, exponent.rdc.v) | 
0000 3 : INPUT PARAMETERS: | 
$000 98 : Both base and exponent are D COMPLEX*16 numbers, each consisting 
0000 (99: of a D REAL*8 real part and a D REAL*8 imaginary part. Both are | 
4 199 : CALL BY VALUE. 
0000 108 : IMPLICIT INPUTS: | 
0909 104 3 NONE | 
0009 106 : OUTPUT PARAMETERS: 
000 107 ; 
0000 109; — 
0000 110 ; IMPLICIT OUTPUTS: 
000 111; 
$000 118-3 — 
9000 116 > FUNCTION VALUE: 
$00 118 : The D COMPLEX*16 (REAL*8, REAL*8) result of taking the 
33 117 ; COMPLEX base to the COMPLEX exponent power is returned 
00 118 ; in registers RO-R3. This is a violation of the VAX 
4 119 ; calling standard, but is excused for compiled code 
: y 3 support routines. 
9 j > SIDE EFFECTS: 
000 ! : ; Modifies registers RO-R3. 
oS ' § : Possible error signals are: 
99 1 3 : MTHS_INVARGMAT if base is (0.,0.). 
0 129 ; MTH$S_FLOOVEMAT if floating overflow occurs. | 
0 1350 ; MTHS_SINCOSSIG if absolute value of the imaginary part of 
131; (exponent * CLOG(base)) > Pl#2*#30. 
: § : SS$_ROPRAND if reserved floating operand is fetched. 


* 
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sENTRY OTSSPOWCDCD_R3, “M<> ; disable integer ovflo 
MTHSFLAG_JACKET ; establish math error handler 


6D O0000000'GF 9 MOVAB G*MTHSSJACKET_HND, (FP) 
set handler address to jacket 


handler 


> 


Get complex logarithm of base 


Ss 
S 
3 3 3 2 3 I OOOO OOOO OGOODOOO 


9 
ao 
3 139 
3 jes 
SE 10 (C2 9 146 SUBL2 #16, SP 3 return complex on stack 
AC OF C 14 PUSHAL base(AP) ; address of base 
A SF OO0F 144 PUSHAB 4(SP) ; address of result 
00000000 ' GF 0 FB § 146 CALLS #2, G*MTHSCDLOG : (SP) gets LOG(base) 
9 167 ;¢4 
019 198 : CLOG(base) is at (SP). Multiply by exponent. 
; 4 133 ; Do multiplication out of Line. 
7E iC aC 7D 0019 «151° MOVQ exp+8(AP), =(SP) ; Put exponent on stack 
7E 14 AC 7D 001D 13 MOVQ exp(AP), =(SP ; CLOG(base) is already there! 
00000000'GF 08 FB 0021 15 CALLS #8, G*OTSSMULCD_R3 : RO-R3 gets CLOG(base) * exp 
058 188 ss 
028. 138 : Now compute CEXP(product) 
$58 188 
7E 3 7d 00 139 mova Ree -(SP) 3 put product on stack 
7E 5 7D 0028 160 MOVQ RO, -(SP) 
SE 10 C2 OO2E 161 SUBL wié, SP ; Make room for result 
10 AE 9F 0031 166 PUSHAB 16(SP) 3; Address of product 
04 Ag oF + 4 16 PUSHAB 4(SP} 3; Address of result 
00000000 ' GF 0 FB 0037 164 CALLS #2, G*MTHSCDEXP :; Result is at (SP) 
50 8€ 7D OO3E 165 MOVG (SP)+, RO : Pop result into RO-R3 
52 BE 7D 49d 196 mOVa (SP)+, R2 
04 0044 198 RET z; all done, exit 
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OTSSPOWCOCD a i PS 9 91:33:99 AX/VMS Macro V04-00 Page 5 
Symbol table 6-SEP-1984 11:27:47 (CMTHRTL.SRCJOTSPOWCDC.MAR; 1 (4) 
BASE = 

EXP = 14 

MTHSSJACKET_HND eeereree x 1 

MTHSCDEXP eeentere a 

MTHSCOL geererese x 

OTSSMUL geeeeree ¥ 


CD _R3 
OTSSPOWCDED_R3 00000000 RG 1 


' 
+ eeeeoeoeoeeoooeaeecon oo 


PSECT name Allocation PSECT No. Attributes 
- ABS, 00000000 ( 0.) 9 ( 0.) NOPIC USR CON ABS - LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
-OTSSCODE 00000045 ( 69.) O1¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 


$wowocecccoccesccecccccon} 


Phase Page faults 


Initialization 30 00300: 00- 18 0:00:00.78 
Command processing 11 sb Shs 4 0:00:03.32 
Pass 7 89:82: 0.59 0:00:02.95 
Symbol table sort 0 0:00:00.00 00:00:00.01 
Pass 2 44 00:00:00.48 00:00:02.02 
Symbol table output 00:00:00.01 es ea 
Psect synopsis output 00:00:00.03 0:00:00.05 
Cross-reference output § es $8 4 00:00:00.00 
Assembler run totals 27 00:00:01.96 00:00:09.15 


The working yee Limit was 750 pages. 
tt bytes {3 pages) of virtual memory were used to buffer the intermediate code. 
ere were 


pages of symbol table space allocated to hold 7 non-local and 0 local symbols. 
229 source Lines were read in Pass 1, producing 11 object records in Pass 2. 
1 page of virtual memory was used to define 1 macro. 


+ 
; Macro Library statistics H 


Macro Library name Macros defined 


-$255$DUA28: CSYSLIBISTARLET.MLB; 2 0 
0 GETS were required to define 0 macros. 


| 

| 

| 

Sipceticcaaaag ; 
Psect synopsis ! 


There were no errors, warnings or information messages. 
MACRO/ENABLE =SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:O0TSPOWCDC/OBJ=OBJ$:OTSPOWCDC MSRC$:MTHJACKET/UPDATE=(ENHS :MTHJACKET) MS 
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